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Abstract 

We analyze the observational and theoretical constraints on "Einstein-aether theory" , a generally 
covariant theory of gravity coupled to a dynamical, unit, timelike vector field that breaks local 
Lorentz symmetry. The results of a computation of the remaining post-Newtonian parameters are 
reported. These are combined with other results to determine the joint post-Newtonian, vacuum- 
Cerenkov, nucleosynthesis, stability, and positive-energy constraints. All of these constraints are 
satisfied by parameters in a large two-dimensional region in the four-dimensional parameter space 
defining the theory. 
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Alternative theories of gravity that deviate from general relativityhave been ruled out or 
severely constrained systematically as observations have improved y], Q . At this stage the 
most studied surviving alternative is "scalar-tensor theory", of which Jordan-Brans-Dicke 
is a well-known example. This sort of theory is a simple extension of general relativity 
containing a fundamental scalar field. The existence of such scalar fields has been suggested 
by moduli fields that arise in higher dimensional approaches to gravity and quantum gravity, 
such as string theory. 

Vector-tensor theories consisting of general relativity coupled to a dynamical vector field 
are much more complicated than scalar-tensor theories, due to the metric derivatives that 
appear in the covariant derivative of the vector field. They also appear to suffer from a glaring 
problem: because of the indefinite signature of the spacetime metric, some of the degrees 
of freedom are always associated with negative energies. This problem need not occur, 
however, if the vector field is constrained to have a fixed timelike magnitude. Such a vector 
field specifies a particular rest frame at each point of spacetime; hence, it "spontaneously" 
breaks the local Lorentz symmetry in a dynamical fashion. These theories — of which there 
is a four-parameter family — preserve the full diffeomorphism symmetry group of general 
relativity, and gravity is still described by the curvature of the spacetime metric. The 
existence of such a vector field is not as well-motivated as a scalar field, however a number 
of approaches to quantum gravity have very tentatively suggested that Lorentz symmetry 
might be broken. If general covariance is to be preserved, then any Lorentz violating vector 
(or tensor) field must be dynamical. 

In this paper we discuss the complete collection of currently available observational and 
theoretical constraints on unit-vector-tensor theories, combining our new result for the re- 
maining post-Newtonian parameters with previously established constraints. Surprisingly, 
all of these constraints are compatible with ranges of order unity for two coefficients in the 
Lagrangian. We are aware of no other theory that comes this close to so many predictions 
of general relativity and yet is fundamentally different. 

We call the timelike unit vector the aether, and the coupled theory Einstein-aether theory, 
or ae-theory for short. A review of its properties together with a collection of references to 
earlier work was given in Ref. jj]- Ae-theory is defined by the action for the metric g a b 
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and aether u a . The most general action that is diffeomorphism-invariant and quadratic in 
derivatives is 

S = f d A x^(R + K ab mn V a u m V b u n + \{u a u a - 1)) (1) 



where 



K ab mn = ci g ab g mn + c 2 5 a J b n + c 3 + c 4 u a u b g mn . (2) 



The coefficients 01,2,3,4 are dimensionless constants, R is the Ricci scalar, and A is a Lagrange 
multiplier that enforces the unit constraint. The expression R ab u a u b is proportional to the 
difference of the C2 and C3 terms via integration by parts, hence is not an independent term. 
Note that since the covariant derivative of u a involves the Levi-Civita connection, which 
involves first derivatives of the metric, the aether part of the action in effect contributes also 

to the metric kinetic terms. We adopt the metric signature (H ), and units are chosen 

such that the speed of light defined by the metric g ab is unity. Other than the signature 
choice we use the conventions of Ref. |4(. 

Observations have already severely constrained Lorentz symmetry violation in the matter 
sector 0, E|, hence to a very good approximation matter must couple universally to one 
metric, which we take to be g a b. Our goal is to determine the observational and theoretical 
constraints on the dimensionless parameters q. 

In the weak-field, slow-motion limit ae-theory reduces to Newtonian gravity j2j, with a 
value of Newton's constant Gn related to the constant G in by 

G v = G(l-^)~\ (3) 
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Here we have written cu for c\ + C4, a notation we generalize below to C123 = c\ + c 2 + C3, 
etc. 

A standard way of beginning to compare an alternative gravity theory to general relativity 
is to examine the first post-Newtonian corrections. For a general metric theory of gravity 
there are ten 'parametrized post-Newtonian' (PPN) parameters characterizing the 

lowest order effects in v 2 /c 2 and dimensionless gravitational potential G^M/c 2 r. Five of 
these parameters, (1,(2, C3; C4; an d 03, vanish identically for any 'semi-conservative' theory, 
i.e. one derived, like ae-theory, from a covariant action principle. Two others, known as the 
Eddington-Robertson-Schiff parameters (3 and 7, characterize respectively the nonlinearity 
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and the spatial curvature produced by gravity. It was previously shown in Ref. |8( that in ae- 
theory /3 — 7 = 1, just as in general relativity. Of the remaining three PPN parameters, two, 
«i, «2, characterize preferred frame effects, and the third, £ (sometimes called the Whitehead 
parameter), characterizes a peculiar sort of three-body interaction. The parameter a 2 for 
ae-theory was computed in Ref. |9( to lowest nontrivial order in the parameters q. 

Here we report on an ab initio computation of all the PPN parameters that confirms 
the previous results and determines the exact values of a.2 and the previously undetermined 
parameters <X\ and £. The parameters are defined by a weak-field expansion of the metric 
in terms of a set of generalized potentials defined by integrals over a fluid source. To 
determine the parameters, one solves the approximate field equations with the fluid source 
in a standard coordinate gauge. The computation is straightforward but lengthy, so the 
details are relegated to an appendix. 

Our results indicate that the "time-time" and "space-space" components of the metric 
are the same in ae-theory and GR to calculated post-Newtonian order, where we refer to a 
nearly globally Lorentz coordinate system specialized to the standard PPN gauge and with 
the aether aligned with the time direction at zeroth-order. The "time-space" components 
of the metric g 0i , i = 1, 2, 3, differ as 



{goi) ae - {goi)GR = — 2 — Vi + Y 



(4) 



where aip are the PPN parameters given explicitly below, while the components of the 
aether are 



u 



u° = l + U 

2 Cl + 3c 2 + c 3 + * {Vi _ ^ + (2-^-03) ^ + m) 



2c 



L23 



2,C\ — -— c$ 



The potentials U, Vi, and Wi are defined by 



(5) 
(6) 



K(t,x) 



} =G N d 6 y 



73„. P 



l x ~ y 



V 

x -y| 2 



(7) 



where p(t, y) is the rest-mass energy density of the fluid source, v y) is the fluid four- 
velocity, and the repeated indices are summed over. 
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The components of the perturbed metric show that the ae-theory PPN parameters are 
given by 



7 = /3 = l (8) 
f = Cl = C2 = C3 = C4 = « 3 = o (9) 



ai = ; 8(ci+ 2 ciC4 i (10) 

2ci -cf + 4 
(2ci 3 - c M ) 2 



«2 - ,„ v 

Cl23(2 - Cu) 

12c 3 c 13 + 2c lCl4 (l - 2c 14 ) + (c? - cl)(4 - 6c 13 + 7c u ) . ( ( ) 



(2-ci4)(2ci-c? + c§) 

Note that the parameters oti and «2 are both of linear order in q when the coefficients are 
small compared to unity and the ratios amongst them are of order unity. 

It is evident from the form of the metric and aether perturbations that the cases C123 = 0, 
C14 = 2, and 2ci — cl + c| = are special, since a% and/or a 2 diverges. Presumably the 
post-Newtonian approximation is not valid when the coefficients are close to these values. 
From the wave speeds (j!5|) below we see that the spin-0 speed vanishes in either of the first 
two cases and the spin-1 speed vanishes in the last case. This corresponds to the absence of 



spatial gradient terms in the action l9j. The case C123 = corresponds j8j to the vector-tensor 
theory of Hellings and Nordtvedt [lO| if the unit constraint on the aether is dropped. This 
theory was shown by Will 1] to be dynamically over-determined and hence observationally 
unacceptable. 

The current best constraints on the preferred frame PPN parameters are |«i| < 10 4 
(from an orbital polarization effect bounded by lunar laser ranging and binary pulsar ob- 
servations) and \a 2 \ <4x 10~ 7 (from a spin precession effect bounded by the alignment of 
the solar spin with the ecliptic). These two conditions can be met with two unrestricted 
parameters to spare, since we begin with four free parameters 0^2,3,4- At the lowest resolu- 
tion, we can just impose a.\ = a 2 = 0. The condition a% = implies c 4 = —c\jc\. (The 
case Ci = can be excluded since it yields cui = —8, which lies far outside the observational 
bound.) Having put ol\ = in this way, it can be seen (the seeing made easier with the help 
of Mathematica) that a 2 can be put to zero in two ways. One way is with C13 = 0, which 
together with the previous condition implies also Cu = 0. This case is degenerate, and is 
briefly discussed at the end of this paper. The other way to put both a\ and a 2 to zero is 
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to determine c 2 and C4 in terms of c\ and C3 via 

c 2 = (-2c? - C1C3 + dj)/3ci 

2/ (12) 

C 4 = -C3/C1 

Thus t/iere is a two-parameter family of ae-theory Lagrangians for which all the PPN pa- 
rameters are identical to those of GR. 

We now consider the other constraints on ae-theory. In alternate gravity theories includ- 
ing Brans-Dicke theory, the Newton constant G*n need not be constant in time. Obser- 
vational bounds on G/G then constrain the theory. In the case of ae-theory, we have the 
relation (J3J), hence Gn is always constant. 

Another constraint arises from the possible discrepancy between Newton's constant and 
the gravitational constant occurring in the equation for the dynamics of the cosmological 
scale factor. In general relativity, the scale factor satisfies the Friedman equation, which 
involves Newton's constant. In ae-theory, when the metric has the standard cosmological 
form (Robertson- Walker symmetry) and the aether is aligned with the cosmological rest 
frame, the aether stress tensor can be constructed purely from the spacetime metric with 
two derivatives, and must be identically divergence free. It must therefore be a linear 



combination of the Einstein tensor G„h and a tensor constructed with the spatial curvature 

e flQ: 

--aether C 13 + 3c 2 



1 Wh a- 

scalar ^R, which turns out to be |7I. Illf: 



T, 



ab 2 



G ab -\^R{g ab + 2u a u b )\. (13) 
J 



The effect of the cosmological aether is thus to renormalize the gravitational constant and 
to add a stress tensor of perfect fluid type that in effect renormalizes the spatial curvature 
contribution to the field equations. The renormalized, cosmological gravitational constant 
is given by [3] 

G cosmo = G(l + C -^±^Y\ (14) 
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Since this is not the same as Gn ©, the expansion rate of the universe differs from what 
would have been expected in GR with the same matter content. The ratio is constrained 
by the observed primordial 4 He abundance to satisfy \G cosrao /G^ — 1| < 1/8, which imposes 
a constraint on the constants q |7(. Remarkably, if the constants are restricted by (|12j) so 
that 0:1,2 = 0, then Gn = G COS mo- Primordial nucleosynthesis then imposes no additional 
constraint. 
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Even when the two gravitational constants coincide, the "curvature fluid" term in (|13|) 
represents a deviation from the Friedman equation in GR if the universe has non-zero spatial 
curvature. Observations have shown that the spatial curvature must be very small today, 
and it would have been even less important in the past when the relative contribution 
of matter and radiation would have been even more important. It thus seems unlikely 
that an interesting constraint can be obtained from this term. Another potential source of 
cosmological constraint is the modification of the primordial fluctuation spectrum [12], but 
this has not yet been worked out in full detail. 

A further constraint on ae-theory comes from the possibility that the gravity and aether 
waves travel at less than the speed of "light" — that is, less than the limiting speed de- 
termined by the metric governing the propagation of matter fields. In this case, high 
energy matter moving inertially through the vacuum would produce vacuum Cerenkov ra- 
diation of gravitational and aether shock waves. A detailed analysis of this process and the 
corresponding observational constraints from ultra-high-energy cosmic ray observations was 
carried out in Ref. j^]. The constraints are characterized by very small numbers, ranging 
between 10~ 15 and 10 -31 , depending on the wave-mode type and emission process. These 
are all one-sided constraints, since they apply only when the wave speeds are smaller than 
the speed of light. To a first approximation then, the constraints imply that the wave speeds 
must be greater than the speed of light. 



Some authors 



12L have suggested that superluminal propagation be excluded a 
priori on the grounds that ae-theory should be viewed as an effective description of an un- 
derlying Lorentz invariant theory in a configuration with broken Lorentz symmetry. How- 
ever, this is a purely theoretical bias, with no observational basis that we can see. Moreover, 
superluminal propagation does not threaten causality, as long as there is a limiting speed 
in at least one given reference frame, as there is in ae-theory. We thus adopt a phenomeno- 
logical stance, allowing for superluminal propagation unless-and-until it is observationally 
ruled out. 

There are five gravitational and aether wave modes in ae-theory: two correspond to the 
usual gravitational spin-2 waves, two are a transverse spin-1 aether-gravity wave, and one is 
a longitudinal spin-0 aether-gravity wave. The squared speeds of these modes are determined 
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by the constants q, and are given by 



spin-2 1/(1 -ci 3 ) 

spin-1 (ci - -c? + -Cg)/ci 4 (l - C13) (15) 

spin-0 ci 23 (2 - ci 4 )/ci 4 (l - ci 3 )(2 + C13 + 3c 2 ). 

If we impose the 0^,2 = conditions (|12jl . the Cerenkov constraint that the spin-2 and spin-0 
wave speeds be superluminal restricts c\ and c 3 to the region 

< C13 < 1 

(16) 

< (ci - c 3 ) < ci 3 /3(l - ci 3 ). 
These conditions also ensure that the spin-1 wave speed is superluminal. 

In addition to the observational constraints already mentioned, there are two theoretical 
constraints, coming from the requirement that the wave modes be stable — i.e. have real 
frequencies — and that the energy of the modes be positive. The first requirement is already 
guaranteed by the no-Cerenkov condition that the speeds be greater than unity. The signs 
of the energy densities of the wave modes, averaged over a cycle, are given by 15| 

spin-2 1 

spin-1 (2c a -c? + cg)/(l-c 13 ) (17) 
spin-0 Ci 4 (2 — C14). 

The spin-2 modes always carry positive energy. If the <x\p, = conditions (fT2j) and the 
superluminal conditions ()16j) are satisfied, then we find that the spin-1 and spin-0 modes 
also carry positive energy. By contrast, if the speeds are sub-luminal, then the latter two 
modes carry negative energy. Thus, not only the Cerenkov constraints, but also the positive 
energy requirement excludes the case of sub-luminal wave speeds. 

We earlier pointed out that an alternate way to set a.\ = a 2 = is if Ci 3 = Ci 4 = 0. In 
this case G^/G cosmo = (1 + 3c 2 /2), so nucleosynthesis would impose a constraint on c 2 . A 
more serious concern is that the spin-0 and spin-1 wave speeds (115)1 diverge in this case, 



because there are no time derivative terms in the aether field equation [1J|. At the same 
time the energy density ([T7Jl of the spin-0 mode goes to zero, but that of the spin-1 mode 
remains finite. We have not worked out the observational signatures of this behavior. 

It is rather non-trivial that the PPN parameters are identical to those of GR and that the 
vacuum Cerenkov, nucleosynthesis, stability, and positive energy constraints are all satisfied 
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in a large two-dimensional region (jl2H6|) in the four-dimensional q parameter space. To 
further constrain the parameters one should look to strong-field effects or radiative processes. 
Since the gravitational radiation damping of binary pulsars has been found to agree with 
GR to one part in 10 3 , one would expect constraints of order 10~ 3 from such observations. 
The ae-theory radiation damping rate has been determined to lowest non-trivial order 
but higher-order terms are still needed to set definitive constraints. 

One particular strong-field effect that can be examined is the existence and nature of 
black hole solutions to the vacuum field equations. Some alternate theories of gravity whose 
PPN parameters are equal or close to those of GR do not admit regular black hole solu- 
tions. In these theories, astrophysical collapse would produce something other than a black 
hole — perhaps a naked singul arity , or a bounce — which may not be difficult to rule out ob- 



servationally. Recent studies [12 1 have shown, however, that ae-theory does admit regular 
black hole solutions. 

In conclusion, ae-theory provides a healthy sparring partner for GR, a role that was 
previously played only by scalar-tensor theory. One of the motivations for exploring this 
theory is the idea that the existence of a preferred timelike direction in spacetime could play 
a role in the solution to fundamental problems such as the nature of dark energy or quantum 
gravity. Although it is likely that radiation and strong-field constraints will limit all of the 
theory's parameters to be small compared to unity, they would still allow the aether field to 
exist and possibly play such a role. 
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APPENDIX: CALCULATION OF AE-THEORY PPN PARAMETERS 

This appendix provides details of the calculation of the Parameterized Post-Newtonian 
(PPN) parameters, ai,a2,a3,/3,7,C:L,C2,C3>C4>£> f° r ae-theory. The PPN formalism is defined 
in a weak-field, slow-motion limit, and describes the next-to- Newtonian order gravitational 
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effects in terms of a standardized set of potentials and these ten parameters. We will 
determine the PPN parameters by solving the field equations with a perfect-fluid source 
in a standard coordinate gauge order-by-post-Newtonian-order. More detailed explanations 
of theprocedure and the general PPN formalism can be found in the classic reference of 
Will |J. 

We will follow the conventions of [4], with the following exceptions and additions. The 
metric will have signature (+, — , — , — ). We will write the approximate equations in terms of 
components of tensors with all indices lowered. Spatial indices will be indicated by lowercase 
Latin letters from the middle of the alphabet: i,j,k,.... Spatial indices are raised and 
lowered with the spacetime metric; e.g. t>j = v a g a i = —v l + (higher order terms). Repeated 
spatial indices are summed over. We will write, e.g. c u for c\ + c 4 , etc. Note also that we 
follow modern conventions for denoting the post-Newtonian order of quantities (see, e.g. 
which differ from the conventions of (l | . 

The ae-theory field equations follow from the action (JI|) . with an additional perfect-fluid 
source coupled in the standard way to the metric g ab , and uncoupled to the aether u a . We 
have the Einstein equations, written in the non-standard form 

R ab = (S cd + MGT cd ) (5 c a 5 d b - ^g ab g cd ), (A.1) 



where 



Sab = Vm(j( a m Ub) — J m (a U b) ~ J(ab)U m ) 

+ ci (V m Ua V m u b - V a u m V b u m ) + c A u a ii h (A. 2) 



with 



and 



also 



+ \u a U b + -gab{J n rr ynU % 



J a m = (cig ab g mn + c 2 5 a J h n + c 3 5 a J b m + c iU a u b g mn )V b u n , (A.3) 
if = u b V b u a ; (A.4) 



T ab = (p + P U + p)v a v b - pg ab , (A.5) 

where v a is the four-velocity, p the rest-mass-energy density, IT the internal energy density, 
and p the isotropic pressure of the fluid. We also have the aether field equation 

V a J a m - c 4 u a V m u a = Xu m ; (A.6) 
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and the constraint 

g ab u a u b = 1. (A.7) 
Eqn. ()A.6|) can be used to eliminate A, giving 

A = u m V a J a m - c 4 u a ii a . (A.8) 

We assume a nearly-globally-Lorentzian coordinate system and basis with respect to 
which, at zeroth order, the metric is the Minkowski metric r) a b and the aether is purely 
timelike. The fluid variables are assigned orders of p ~ II ~ p/p ~ (v 1 ) 2 ~ 0(1). Taking 
the time-derivative of a quantity will effectively raise its order by one-half: X ~ 0(N) — > 
dX/dt ~ 0(N + 1/2). We assume that the components of the metric perturbations h a b with 
respect to this basis will be of orders 



'mi 



0(1) + 0(2), ^-0(1), hoi~ 0(1.5). (A.9) 



This assignment preserves the Newtonian limit while allowing one to determine just the first 
post-Newtonian corrections. The aether perturbations 5u a are assumed to be of orders 

5u° ~ 0(1), 5u l ~ 0(1.5). (A.10) 

Lower orders are disallowed by the field equations, given the above orders of h a b- We will 
assume that h a b and Su a satisfy boundary conditions such that they vanish at spatial infinity. 

The metric components are to be expanded in terms of particular potential functions, 
thus defining the PPN parameters: 

<7oo = 1 - 2U + 2(3U 2 + 2^ w - (2 7 + 2 + a 3 + Ci - 2£)$i 

- 2(3 7 - 2(3 + 1 + C 2 + 0$2 - 2(1 + Cs)<f 3 - 2(3 7 + 3( 4 - 2^)$ 4 + (Ci - 2CM, 



0y = -(l + 2 7 LO* 

^ = ^(4 7 + 3 + en - a 2 + Ci - 20^ + -(1 + « 2 + Ci + 20Wi. 

(A.ll) 

The potentials are all of the form 

F(x) = G N / d 3 yp^-, (A.12) 
J \x-y\ 
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where Gn is the current value of Newton's constant, which we determine below in terms of 
G and the Cj. The correspondences F : f are given by 



U : 1 $1 : ViVi $ 2 : U $ 3 : II $ 4 : p/p 

^: /^p(,) | X ~^ f ( ,^^ - ( "-^ ) (A.13) 
|x — V |x — z\ \y — z\ J \x — y\ 2 



V, : v l Wi 



,J j( x j -yj)(x % -y l ) 



\x-y 



2 



Note that for U, $1,2,3,4, and Vi, 

F >ii = -^irG N pf. (A. 14) 

We will also make use of the 'superpotential' x- 

X = -G N J d 3 y p\x - y\, (A.15) 

which satisfies 



We also note the relation 



X,u = -2U. (A.16) 



X,io = Vi-Wi, (A.17) 

which follows from the formula 

^ J d 3 yp(y,t)f(x,y) = J d 3 yp(y, ty(y, t)^[l + O(l)], (A.18) 
which follows from the continuity equation for the fluid 

P,o+(pv i ),i = 0, (A.19) 

assumed to hold to 0(1.5). 

n 

These potentials satisfy certain criteria of "reasonableness" and simplicity (see [1|, Sec. 
(4.1) for details), and are general enough to describe all known viable theories of gravity. In 
particular, they suffice for ae-theory. The criteria permit goo to depend also on the potential 
X oo, an d gij to depend on x,ij- Such terms, however, can always be eliminated [l| by a 
suitable coordinate transformation that preserves the zeroth-order form of the components. 
The 'standard PPN gauge' is thus defined as that post-Newtonian coordinate frame in which 
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all dependence on x j0 o and x,ij has been removed from, respectively, goo and gij. This fixing 
determines the coordinate frame up to necessary order so that the standard form of the 
metric components is unambiguous. 

In carrying out the calculations, we shall impose the following gauge conditions: 

hijj = -^{h 00 ,i - h jjti ) (A.20) 
h 0i ,i = -3l7,o +0n i4 (A.21) 

where 9 is an arbitrary parameter and ni = U{ — h^. These conditions are suggested by the 
standard conditions for general relativity. As we shall see, the conditions (|A.20|) suffice to 
put gij in standard form, while the fourth condition (jA.21j) standardizes #oo when 

(ci + 2c 3 - c 4 ) , , 

6 ~ (2-c 14 ) • (A ' 22) 

The solving procedure is as follows: 

Step 1: Solve the constraint (jA.7j) for u° to 0(1); 

Step 2: Solve the "time-time" component of the Einstein equation (jA.l|) for g 00 
to 0(1); 

Step 3: Solve the "space-space" components of (jA.l|) for to 0(1); 
Step 4: Solve the "space" components of the aether field equation (jA.6|) for u % 
to 0(1.5); 

Step 5: Solve the "time-space" components of (|A.1J) for g 0i to 0(1.5); 
Step 6: Solve the "time-time" component of (jA.ljl for g o to 0(2). 

The cases in which ci 23 = 0, c u = 2, or 2ci — c\ + c\ = are special in that the found 
solutions diverge. We will presume that the post-Newtonian approximation is not valid in 
these cases, and assume below that they do not hold. See the main text for more discussion 
of this point. 

1. u° to 0(1) 

Solving the constraint (|A.7|) gives 

u° = 1- (l/2)/i 00 (A.23) 
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to 0(1). For the components of u a , we have 

u Q = u a g a0 = 1 + -hoo, (A.24) 

and 

Ui = u a 9ai = -u l + hi = rii + h 0i . (A. 25) 

For later convenience, we will now express the covariant derivatives of u a . The con- 
straint (|A.7jl implies that 

V a u = (A.26) 

to 0(2). Also to 0(2), we have 

V Mj = ~hao,i{l ~ 7^00) + h m ,o + n ifi , (A.27) 



and 



Ui = u°V Ui = ~hao ti (l - h 00 ) + hoi,o + n i)0 . (A.28) 



To 0(1.5), we have 

VjUi = riij + -hyfl + h 0[iJ] . (A.29) 



2. 5 oo to 0(1) 

We now solve the "time-time" component of the Einstein equation (|A.1|) for goo to 0(1). 
For the components of R 00 , we have 

RoO — 2^00,ii + ~ {hiO,i ~ 2^**>o) o — ^00,i^00,i + T^QOJ (2/ii7,j — huj) (A. 30) 

to 0(2). At 0(1), we have 

-Roo — 2^oo,m, (A.31) 

Too = p, Tij = 0, (A.32) 

SoO = Jo™,™ ~~ ^00,0 = —Joi,i = ~C\a{^ QUi) ,i = ^00,ii; (A. 33) 

S i:i = 0. (A.34) 

The field equation becomes 

(1 - ^) V« = 87rG >> (A. 35) 
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which gives /loo to 0(1), 

h 00 = -2U, (A.36) 



with Newton's constant 



G N =(l- C -fY l G. (A.37) 



3. gtj to 0(1) 

We now solve the "space-space" components of ()A.1|) for to 0(1). We have to 0(1) 

Rij ~zhij,kk ~Z^kk,ij ^fc(ij')fc ~zh()0.ij 

I (A.38) 

where we have imposed the gauge condition f)A.20|) in the second step. Using (|A.32|) . ()A.33j) . 
and (|A.34J) . the field equation becomes 

h i5 ,kk = '&'KG N p5 i j, (A.39) 

giving 

tHj = -2US i:j . (A.40) 



4. u* to 0(1.5) 

We now solve the "space" components of the aether field equation (jA.6|) for u l to 0(1.5), 
making use of the gauge condition (|A.20|) and our earlier results (jA.24j) , (|A.36|) , and (|A.40|) . 
At 0(1.5) equation (jA.6|) has the form 



To 0(1.5), 



and 



J\a = Joi,o - J ji;j = 0. (A.41) 



Joi,o — ci4(Voitj) i o — — —hoo,io — — ~^~Xfiijj) (A. 42) 



= CiTlijj + C2J,rijji + - (2c_/lQ[ij]j + (c + + 3C2)x,0ijj) : 



(A.43) 
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where c_ = c\ — C3, c + = C13. The aether field equation can then be written 
/ c i \ c 

(can* + y/ioi + -ji^ c i + 3c 2 + c 3 + c 4 )x,ioJ . . - [-^hajj - 023^,3)^ = °- (A.44) 
Taking the spatial divergence of the left-hand side gives the relation 

n^ijj = Axflujj, (A.45) 

where 

A = - 2Cl + 3C2 + C3 + C4 , (A.46) 

2Ci23 

which we can solve for n^. Substituting into ()A.44|) . imposing the gauge condition (|A.21j) . 
and using our earlier results, we solve (|A.44|) : 

m = -u l = - J-( c _/i 0l - (2c x A + c_(^ + Ae))x >Qi ). (A.47) 
5. g 0l to 0(1.5) 

We now solve the "time-space" components of (jA.ljl for to 0(1.5), making use of the 
gauge conditions (|A~20j) and (TOT|) and the earlier results (TOIJ) . (TOB) . (TOn|l . and (|A~47jl . 

We have to 0(1.5) 



Rot — ho[i,j]j + hj\j t i]Q 

= \(h Q i+ l -(l-2Ae) X> ^j P 



(A.48) 



Also to 0(1.5), 



and 



We have 



T 0i = pvi, (A.49) 



1 1 

Soi = —J(io),o + 2^ m ,m. = — 2 (^ 0i >° ^.i)" (A. 50) 



U,3 



C 3 n i,j),j + g (2c_/io[j,i]jr' + c +hij,jO + C2hjj t io) (A. 51) 



where 



E = (c\ + 3c^ + 4cic 4 - 2c_c + Ad^j . (A.52) 
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With (lA~42j) . this gives 

Soi = (^^hoi + —x,oi) • (A.53) 

v 4Ci I / ,jj 

The field equation becomes 

(1 - V# = WnGpoi + (E + A9- -)x,oijj, (A.54) 



giving 



_n c - c +\-i 

oi — U — 



2ci 

6. 5oo to 0(2) 



({E + Ae- l -) X & + ^- C -f)v). (A.55) 



We now solve the "time-time" component of (jA.l)) for goo to 0(2), making use of the 
gauge conditions (fA~20jl and (fA~2T]l and the earlier results (fA~24|) . fA~3fijl . (fA~40jl . (fA~47j) . 
and ()A.55|) . Define /i o — #00 — 1 + 2Z7. Then, from eqn. ()A.30|) . we have 

Roo = \ (hoo -2U- 2U 2 + 8$ 2 - 2A6 X , 00 ) ^ (A.56) 



Also, 



T 00 = p(l + n + ^-2f/), (A.57) 
= pviVj + pSij. (A. 58) 

9m{T ab g ah ) = (1 - 2f/)(T 00 (l - 2Z7) - r i (l + 210), (A.59) 



so that 



Too - ^oo( W) = + n + 2{v iVi - U)) + |p 

' y -(f/ + 2$ 1 -2$ 2 + $3 + 3$ 4 ) 4r 



(A.60) 



8nG 

To aid the reader in sorting through the terms appearing in Sab, we note that to 0(2), in 
the chosen gauge, 

3 

V lUl = (- + A)x,oii (A.61) 

and 

(V «i),i = -\ ( - 2U + /ioo - U 2 - 2 + (1 + 0)A) X ,oo) ■ (A.62) 
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After some bookkeeping, we find that 



and 



S 00 = ^ ( - 2U + h 00 - -U 2 + 9$ 2 ) - c 14 (- + (1 + e)A)x,wu, (A.63) 



Su = °-f{\u 2 - $ 2 ),« - (c+ + 3c 2 )(^ + A) Xmt . (A.64) 



We thus have 



Soo ~ 29ooSab9 ab — 2 (^oo + Su) + 2U S^, 



C14 

T 
1 /,3 



■ f - 2Z7 + /ioo - 2f/ 2 + 8$ 2 ) (A.65) 
- ((- + A){2c x + 3c 2 + c 3 + c 4 ) + c 14 A^x,oou- 
We combine 1A.56|) . (jA.60|) . and fjA.65jl . solve the field equation, and obtain 

h m = 2U 2 - 4$i - 4$ 2 - 2$ 3 - 6$ 4 + Qx,oo, (A.66) 

where 

Q = (1 - ^i)- 1 ((2 - c 14 )0 + (d + 2c 3 - c 4 ))A (A.67) 
Finally, we move into the standard gauge by choosing # so that Q vanishes: 

9-^9 = - Cl + 2C3 ~ C4 . (A.68) 

z — Ci 4 

7. Summary 

We now collect our results, eqns. (lA~36l) . (IA~40|) . (lA~55l) . and (lA~66l) for the metric compo- 
nents, and (jA.23|) and ()A.47|) for the aether, imposing the gauge conditions (jA.20|) and (|A.21|) 
with 9 = 8q, and using the relation (|A.17j) . For the metric, we have 

g 00 = 1 - 2U + 2U 2 - 4$i - 4$ 2 - 2$ 3 - 6$ 4 (A.69) 
^• = -(l + 2f/)<% (A.70) 



2c 



2C i 2 ((E + A6 - \ + 2(2 - c 14 ))^ -(E + A6 - ~)W;) . (A.71) 
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We can extract the PPN parameters by comparison with the standard forms (jA.ll|) . We 
find 



7 = (3 = 1 (A.72) 

f = Ci = C 2 = Cs = C 4 = « 3 = (A.73) 

ai = - o 8(ci+ 2 ClC4 l (A.74) 
2ci - cf + c§ 

(2ci 3 - c i4 ) 2 

«2 ~~ 



ci 23 (2 - cu) 

12c 3 c 13 + 2c lCl4 (l - 2c 14 ) + (c 2 - c|)(4 - 6c 13 + 7c 14 ) 



(2-c 14 )(2 Cl -cf + c|) 



For the aether, we have 



where 



(A.75) 



w° = 1 + U (A.76) 
u* = -(A + B)V i +(A-.B)W i , (A.77) 



2ci + 3c 2 + c 3 + c 4 (2 - c 14 )c_ 

97 ' Vc ^—72- ( A - 78 ) 

ZC123 /c i — c i + c 3 
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